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ABSTRACT 

Hydraulic  ram  is  the  dynamic  loading  of  fuel  tanks  when  impacted  by 
bullets  or  other  projectiles 0  During  impact  and  penetration  of  the  fuel  cell, 
intense  pressure  waves  are  generated  by  the  projectile. 

A  ballistic  range  was  built  and  experimental  testing  was  conducted 
to  study  hydraulic  ram  phenomena.    A    Q.22  caliber  rifle  was  used  to  accel- 
erate projectiles  at  velocities  in  the  range  of  0.38  km/sec  into  a  trans- 
parent, water  filled  tank.    Shape  and  intensity  of  the  shock  wave  pressure 
pulse  induced  as  a  result  of  projectile  impact  were  determined  using  a  dual 

shadowgraph  system.     Peak  pressures  were  found  to  be  as  high  as    ^.50  x 

3  2 

10     kg/cm    .      The  rate  of  energy  transfer  to  the  fluid  by  the  projectile  was 

determined  experimentally  and  compared  with  analytical  predictions.    A 

characteristic  time  was  defined,  establishing  the  separation  point  between 

the  shock  and  cavity  phases  of  hydraulic  ram. 
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I.      INTRODUCTION 

When  enclosed  fluid  containers  are  impacted  by  high   speed  pro- 
jectiles ,   severe  damage  can  occur  as  a  result  of  the  phenomenon  known 
as  hydraulic  ram.     Numerous  experimental  investigations  of  the  hydraulic 
ram  effect  have  been  conducted  including  research  by  Cole,  Stepka, 
Lundstrom,  and  Yurkovich  (Refs  .   1-7).    An  investigation,  presently  being 
conducted  at  the  Naval  Postgraduate  School,  graphically  demonstrates 
this  effect.    The  damage  shown  in  Tigs.    1  and  2  occurred  when  a  plexi- 
glas  test  tank,  simulating  an  aircraft  fuel  cell,  was  impacted  by  a  2.92 
gram,   0.222  caliber  projectile  traveling  at  a  velocity  of  0.889  kilometers 
per  second  (2910  ft/sec). 

In  determining  aircraft  survivability  the  destruction  due  to  hydraulic 
ram  becomes  of  paramount  importance  and,  in  fact,  is  the  major  cause  of 
damage  when  an  aircraft  is  hit  by  small  arms  fire  or  warhead  fragments 
(Ref.   6).     This  is  further  substantiated  by  the  fact  that  fuel  tank  (or  cell) 
components  comprise  the  greatest  portion  of  the  target  area  in  modern 
attack  and  fighter  aircraft.    There  are  numerous  examples  of  catastrophic 
failure  of  fuel  cells  which  were  subjected  to  ballistic  impact  (Ref.  6). 
Other  kill  modes  due  to  hydraulic  ram  include:    fuel  starvation,  explosion 
and  fire,  and  damage  to  aircraft  components  adjacent  to  the  fuel  cell. 
In  space  flight  the  hazard  of  meteoroid  impact  into  liquid  propellant  tanks 
is  another  area   of  concern  because  an  impact  of  sufficient  energy  could 
cause  catastrophic  bursting  of  the  tank  wall. 
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Fig.   1       Hydraulic  Ram  Damage  to  the  Test  Tank. 
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Concern  with  the  aircraft  fuel  tank  vulnerability  problem  did  not 
begin  with  the  recent  air  conflict  in  Southeast  Asia.     Since  World  War  II, 
there  have  been  efforts  to  reduce  hydraulic  ram  damage  to  an  acceptable 
level;  however,  with  few  exceptions,  the  efforts  have  been  "trial  and 
error"  procedures.    The  result  has  been  that  many  schemes  have  been 
tried,  but  none  has  been  found  fully  effective.    While  not  being  satis- 
factory, many  of  these  ideas  have  shown  promise  in  defeating  at  least 
some  of  the  hydraulic  ram  components.     Several  important  fuel  tank  defense 
schemes  were  developed  and  tested  by  North  American  Rockwell  during  the 
late  1960's.    These  concepts  are  published  as  a  survivable  fuel  tank  sys- 
tems handbook  (Ref.   8)  in  which  the  basic  design  features  are  divided  into 
"intrinsic"  and  "parasitic"  defense  concepts.     Intrinsic  defense  concepts 
utilize  specific  aircraft  configuration  considerations  to  reduce  system 
vulnerability.     Examples  would  be:    tank  shapes  with  smooth  contours, 
low  surface-to-volume  ratios ,  few  corners  and  fittings ,  and  fuel  system 
considerations  to  isolate  the  cell.    This  type  of  defense  concept  is  thus 
"intrinsic"  to  the  basic  airframe  design  and  can  be  effectively  applied 
only  when  the  aircraft  is  in  preliminary  design  stages.     Parasitic  defense 
concepts  include  protective  devices  and  systems  added  to  the  aircraft 
primarily  to  enhance  fuel  cell  survivability;  these  concepts  include  ele- 
ments not  otherwise  essential  to  normal  aircraft  function.    Examples  are: 
self-sealing  tanks  and  the  addition  of  reticulated  foam  inside  fuel  cells 
to  reduce  the  hydraulic  ram  effect. 


13 


The  purpose  of  this  study  is  the  investigation  of  the  hydraulic  ram 
phenomenon  itself.     Fuel  tank  response  to  projectile  penetration  cannot 
be  fully  understood  without  a  prior  comprehension  of  the  individual 
hydraulic  ram  components.    This  understanding,  plus  a  knowledge  of 
events  which  occur  during  the  damaging  process,  must  be  obtained  in 
order  to  change  the  tank  arrangement  in  such  a  way  as  to  properly  influ- 
ence the  damage.     Once  fuel  cells  are  protected  against  this  type  of 
threat,  overall  aircraft  vulnerability  will  be  lessened. 


14 


II.     NATURE  OF  THE  PROBLEM 

It  has  been  found  (Refs  .   6  and  7)  that  the  loading  and  deformation 
of  fuel  tanks  by  hydraulic  ram  is  a  combination  of  a  number  of  events. 
These  events,  known  as  the  hydraulic  ram  components,  are  outlined  in 
Figs.  3  and  4.    The  column  on  the  left  in  Fig.  3  describes  the  motion  of 
the  projectile  as  it  passes  through  the  tank.    When  the  projectile  enters 
the  tank,  a  shock  wave  is  formed  and  cracks  may  develop  due  to  the  high 
pressure  and  stresses  near  the  entry  point.    As  the  projectile  traverses 
the  fluid,  high  pressure  is  generated,  and  energy  is  imparted  to  the  fluid 
through  projectile  drag.    This  increased  energy  sets  the  fluid  in  motion 
forming  a  cavity.     Both  the  fluid  motion  and  the  subsequent  cavity  forma- 
tion and  collapse  impose  stresses  on  the  tank  walls.     Finally,  if  the  pro- 
jectile has  sufficient  kinetic  energy,  it  will  exit  the  tank.    As  it  exits 
from  the  tank,  a  local  area  of  compressed  fluid  is  formed,  and  cracks  are 
produced  around  the  exit  hole.    As  shown  in  Fig.   4,  the  wall  loadings 
plus  the  cracks  at  the  entry  and  exit  points  can  result  in  various  types  of 
damage. 

Excluding  crack  formation  due  to  the  actual  entry  and  exit  of  the 
projectile,  hydraulic  ram  is  composed  of  two  main  separable  effects;  a 
shock  phase  and  a  drag  or  cavity  phase,  both  of  which  are  important  in 
causing  damage.     During  the  initial  impact,  the  projectile  impulsively 
accelerates  the  fluid,  thus  generating  an  intense  pressure  field  of  short 
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Fig.  3       Outline  of  Hydraulic  Ram  Phenomenon 
adapted  from  Reference  10. 
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duration.    This  phenomenon  is  characteristic  of  the  shock  phase  of 
hydraulic  ram.    The  contribution  of  the  shock-phase  pressure  pulse 
toward  fuel  cell  failure  is  particularly  important  in  the  case  of  small, 
high  velocity  projectile  impact.    After  passage  of  the  projectile,  pressure 
phenomena  become  increasingly  complex  due  to:     ].)  multiple  sources  of 
the  pressure,    2)  reflections  of  the  pressure  waves,    3)  motion  of  the 
fluid,  and  4)  tank  wall  deflection.    These  events  lead  to  the  formation 
of  a  cavity  and  mark  the  beginning  of  the  cavity  phase.    The  formation  of 
a  large  cavity  behind  the  projectile  is  the  most  striking  phenomenon 
associated  with  ballistic  penetration  of  fluids.    The  pressure  field  pro- 
duced by  the  cavity  formation  and  collapse  is  of  low  amplitude  but  long 
duration  and  has  a  dynamic  behavior  loosely  analogous  to  that  of  a  cavity 
formed  by  an  underwater  explosion. 

The  shock  and  cavity  phases  that  have  just  been  described  have 
been  treated  in  this  analysis  as  two  separate  effects.    A  clearer  under- 
standing of  the  relationship  between  the  shock  and  cavity  phenomena  can 
be  gained  by  considering  the  projectile  energy  transfer  time-history  as 
shown  in  Fig.   5  (Ref.  9).     Two  separate  regions  can  be  distinguished 
which  correspond  to  the  two  main  effects  present  in  hydraulic  ram.    The 
first  region  is  characterized  by  a  high  energy  transfer  from  the  projectile 
to  the  fluid.     It  is  while  the  projectile  is  in  this  region  that  the  shock 
wave  is  formed  and  begins  to  propagate.    The  second  region  is  character- 
ized by  a  lower  and  nearly  constant  rate  of  energy  transfer  to  the  fluid . 
This  second  region,  which  lasts  for  the  remaining   time  that  the 
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Fig.   5     Energy  Transferred  to  the  Fluid  by  a 
Projectile  Versus  Time. 
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projectile  is  traversing  the  tank,  is  known  as  the  drag  or  cavity  phase. 
In  this  phase,  projectile  energy  is  transformed  into  the  kinetic  energy  of 
fluid  motion  which  leads  to  the  formation  of  a  cavity. 

The  work  here  was  done  to  further  the  development  of  analytical 
methods  used  in  understanding  the  complex  hydraulic  ram  phenomenon. 
Once  an  ability  to  predict  the  response  of  fuel  cells  to  ballistic  impact 
exists,  overall  aircraft  survivability  will  be  enhanced. 

A.  SHOCK  PHASE 

Shock  wave  formation  and  propagation  due  to  projectile  penetration 
into  a  fluid-filled  tank  will  result  in  very  high  local  pressures  which  can 
cause  catastrophic  failure  of  the  entry  wall.    This  result  has  been  shown 
experimentally  (Refs  .   2-5),  and  it  provides  the  motivation  behind  the 
study  of  the  shock  phase  of  the  hydraulic  ram  phenomenon. 

The  shock  or  impact  phase  can  itself  be  regarded  as  a  succession 
of  interrelated  events  which  will  be  discussed  in  their  order  of  occurrence. 
When  a  fuel  tank  wall  is  impacted  by  a  high  velocity  projectile,  dynamic 
stresses  are  generated  and  transmitted  through  the  wall.    These  stresses 
are  due  to  the  puncturing  and  cratering  action  from  the  projectile  impact 
and  take  the  form  of  compressive  and    circumferential  bending  stress  in 
the  wall  (Ref.  3).    As  the  projectile  is  penetrating  the  tank  wall,  it 
impulsively  accelerates  the  fluid  directly  behind  the  impact  point.    An 
instant  before  impact     the  fluid  is  at  rest,  and  in  the  next  instant  it  is 
traveling  at  the  velocity  of  the  projectile.    This  sudden  acceleration  of  the 
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fluid  by  the  projectile  during  impact  and  penetration  generates  an  intense 
pressure  field  bounded  by  a  shock  wave.    Tank  wall  stresses  are  induced 
by  this  pressure  field  which,  when  added  to  the  dynamic  stresses  due  to 
cratering,  can  cause  catastrophic  failure  of  the  entry  wall.     The  shock 
wave  pressure  field  expands  in  a  hemispherical  shape   away  from  the 
impact  point  with  a  velocity  greater  than  the  acoustic  speed  of  the  fluid. 
Both  the  velocity  of  the  shock  front  and  the  intensity  of  the  related  pres- 
sure field  vary  directly  with  the  amount  of  energy  transferred  to  the  fluid 
by  the  projectile  during  the  initial  impact  (Refs  .   3  and  4).     The  amount 
of  impact  energy  which  is  transferred  to  the  fluid  is  dependent  on  a  num- 
ber of  parameters  including:      projectile  size,  density,  and  geometric 
shape;  tank  wall  material  and  construction;  and  the  physical  properties 
of  the  fluid  (Ref.   5).     In  general,  for  any  given  time  after  impact,  small 
high-velocity  projectiles  lose  a  larger  part  of  their  impacting  kinetic 
energy  than  more  massive,  low  velocity  projectiles  (Ref.  2).    This 
observation  indicates  that  for  the  same  impact  kinetic  energies,  a  more 
damaging  pressure  pulse  for  wall  fracture  will  result  from  impact  by 
small  high  velocity  projectiles  than  with  more  massive  low  velocity 
projectiles  . 

Although  the  peak  pressure  due  to  projectile  impact  is  high,  the 
pressure  field  is  rapidly  weakened  because  of  the  geometric  expansion 
of  the  shock  wave.    An  experimental  example  of  this  pressure  attenuation 
is  provided  by  Stepka  (Ref.   2).     Impact  with  a  0.9  gram  steel  cylinder  at 
an  impact  velocity  of  4.27  km/sec  resulted  in  pressures  of  4.63 
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giganewtons  per  square  meter  (0.67  x  10     lb/in    )  at  a  distance  of  1 . 9 

centimeters  (3/4  inch)  from  the  impact  point.     However,  the  pressure  at 
the  shock  front  decayed  to  ambient  values  within  13  centimeters  (5 
inches)  from  the  impact  point.    These  results  indicate  that  even  for  im- 
pacts  with  high  velocity  and  high  kinetic  energy  projectiles,  the  side  or 
rear  walls  of  any  large  tank  would  not  be  subjected  to  significant  pres- 
sures due  to  the  shock  phase. 

In  summary,  the  shock  phase  of  hydraulic  ram  is  characterized  by 
an  intense  pressure  field  of  short  duration.    It  is  essentially  a  "point 
source"  energy  release  and  results  in  an  expanding,  hemispherically 
shaped,  pressure  wave  front  that  emanates  from  the  point  of  impact.    The 
kinetic  energy  of  the  projectile  is  the  most  important  parameter  in  deter- 
mining the  speed  of  propagation  and  strength  of  the  shock  front.     The 
shock-phase  pressure  pulse  contributes  to  fuel  cell  failure  mainly  in  the 
case  of  small  high  velocity  projectile  impact. 

B.  CAVITY  PHASE. 

After  the  initial  impact  with  the  tank,  which  is  responsible  for  the 
shock  phase,  the  projectile  continues  to  decelerate  as  it  traverses 
through  the  fluid.     It  is  during  this  time  that  the  cavity  phase  of  the 
hydraulic  ram  phenomenon  is  generated  . 

When  an  object  moves  through  a  fluid,  it  feels  a  resistance  known 
as  "drag.."     Drag  forces  are  due  to  an  uneven  pressure  distribution  and 
viscous  effects  of  the  fluid.     If  the  pressure  distribution  is  integrated  in 
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the  direction  of  motion,  the  resulting  force  is  known  as  "pressure  drag." 
Compared  with  the  drag  duo  to  viscous  effects,  pressure  drag  is  the  pre- 
dominant drag  force  (90%)  at  projectile  velocities  (Ref.    10). 

Drag  is  usually  expressed  in  terms  of  a  nondimensional  "drag 
coefficient,"  which  is  dependent  on  the  density  of  the  fluid,  as  well  as 
on  the  size,  velocity,  and  geometric  shape  of  the  object.    A  typical  drag 
coefficient  for  a  high- velocity  projectile  is  0.  15,  but  this  value  will 
increase  by  a  factor  of  ten  or  more  for  a   tumbled  or  flattened  projectile 
(Ref.    10).     In  the  analytical  study  of  hydraulic  ram,  this  coefficient  can 
be  used  to  account  for  the  much  greater  pressure  generating  capability  of 
deformed  and  tumbled  projectiles. 

Applied  to  the  projectile  surface,  the  drag  forces  generate  a  pres- 
sure gradient. between  the  projectile  and  the  fluid.    This  causes  the  sur- 
rounding fluid  to  move  away  from  the  surface  of  the  projectile.    In  this 
manner,  the  projectile  energy  is  gradually  transformed  into  the  kinetic 
energy  of  fluid  motion.     In  contrast  to  the  shock  phase,  the  fluid  is 
accelerated   gradually  rather  than  impulsively  so  the  peak  pressure  is 
lower.    However,  the  duration  of  the  pressure  pulse  is  much  longer. 
Once  a  sufficient  amount  of  energy  is  transferred  to  the  fluid,  it  will 
form  a  void,  known  as  a  cavity,  around  and  behind  the  projectile.    The 
cavity  which  is  formed  would  be  a  complete  vacuum  except  for  two 
phenomena.    It  is  filled  with  fluid  vapors  from  pressure  vaporization 
effects  and  with  entrained  air  from  the  projectile  entry  and  exit  points. 
After  passage  of  the  projectile,  the  hydraulic  ram  cavity  continues  to 
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expand  until  all  fluid  momentum  has  been  dissipated.    The  momentum  is 
lost  to  the  surrounding  fluid  by  fluid  pressures  which  resist  cavity  growth. 
These  same  pressures  including  pressure  due  to  the  force  of  gravity 
cause  the  cavity  to  collapse.     Collapse  of  the  over-expanded  cavity 
would  be  exactly  opposite  to  cavity  growth  except  for  the  presence  of 
fluid  vapors  and  trapped  air.    These  vapors  and  air  act  as  a  spring  to  make 
complete  collapse  impossible.     From  this  point  on  the  cavity  acts  much 
like  an  explosion  cavity;  the  cavity  oscillations  generate  low  amplitude 
but  long  duration  pressure  pulses.    The  cavity  may  oscillate  through 
several  cycles  until  the  energy  of  the  fluid  is  dissipated.    The  pulse  width 
and  peak  pressure  produced  by  the  cavity  oscillation  depend  on  the  total 
energy  transferred  to  the  fluid,  the  amount  of  gas  trapped  in  the  cavity, 
and  the  shape  of  the  cavity  itself  (Ref.  6).     Pressure  pulses  due  to  the 
cavity  phase  of  hydraulic  ram  are  generally  more  damaging  than  those  of 
the  shock  phase  due  to  their  longer  duration. 

The  rate  of  energy  transfer  to  the  fluid  is  the  major  factor  that 
determines  the  severity  of  hydraulic  ram.     Projectile  kinetic  energy  dis- 
sipation is  proportional  to  the  square  of  the  velocity  decay  as  the  pro- 
jectile traverses  the  fluid.     This  indicates  that  the  projectile  deceleration 
time-history  as  it  passes  through  a  fluid  is  of  great  importance   in  the 
study  of  hydraulic  ram. 
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HI .      ANALYTICAL  APPROACH 
TO  CAVITY  OSCILLATION  AND  PRESSURE  FIELD 

Two  analytical  problems  were  investigated  concerning  the  oscilla- 
tory cavity  phase  of  hydraulic  ram.    The  first  problem  involved  the  deriva- 
tion of  a  method  for  predicting  the  collapse  time  for  a  cavity  when  the 
vapor  pressure  within  was  negligible.    The  assumption  applied  to  this 
problem  was  that  the  cavity  contained  a  vacuum.    The  second  problem 
earned  the  same  basic  analysis  even  further  by  modeling  the  oscillation 
of  a  cavity  with  finite  internal  pressures.    This  cavity  was  found  to  behave 
much  like  a  spring  in  its  oscillatory  pattern.     Isentropic  compression  and 
expansion  were  the  basic  assumptions  used  in  the  second  problem. 

A.  COLLAPSE  OF  A  CAVITY  CONTAINING  A  VACUUM  VERSUS  TIME 

Potential  flow  relationships  for  an  incompressible  fluid  were  used 
in  the  formulation  of  this  problem.    The  cavity  was  assumed  to  be  spheri- 
cal in  shape  with  an  initial  radius   RQ.     R  =  R(t)  describes  the  motion  of 
the  fluid -gas  surface  which  travels  at  the  velocity  u.     From  continuity, 
the  flow  field  is  described  by 


<)  r 

where  r  is   the   instantaneous   cavity  radius.      Potential  flow  theory 
along  with  the   momentum   equation   for  an  incompressible   fluid   was 
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used  to  obtain  the  equation  of  motion  for  a   cavity.    It  was    found 
(Appendix  A)    that: 

K  K  *■  T  ft    - (2) 

where  P  (R,  t)   was   the  pressure  in  the  fluid  at  the  cavity  wall.    Assuming 
that   P  (R,  t)  was  known,  equation  (2)  was  manipulated  into  the  general 
form 

f    J  (3) 

For  the  first  problem  with   R  ■     R0  ,    R0  '-     0    and   P  (R,  0)   =    0  ,    the 
integral  yields  the  time  for  the  cavity  to  reach  a  fraction  of  its  initial 
radius 

xH  JX 

IT-  x3 

v      4 

where  X  =  R  /  RQ  .    This  equation  was  numerically  integrated  by  compu- 
ter.   The  result  describes  the  motion  of  a  spherical  cavity  with  time 
when  the  cavity  contains  a  vacuum  and  starts  at  some  initial  radius  RQ  • 
Solutions  for  two  cavities  with  varying  initial  conditions  are  shown  in 
the  computer  outputs  of  Figs.  6  and  7.    These  plots  give  valuable  insight 
into  the  frequency  response  and  an  order  of  magnitude  for  the  time  to 
cavity  collapse.      Figure   6    gives   a   time   to   collapse   of   0.011  seconds 
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Fig.  6     Cavity  Collapse  Time  from  an 
Initial  Radius    of   0.5  Foot. 
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Fig.  7     Cavity  Collapse  Time  from  an 
Initial  Radius   of    1.0  Toot. 
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2 
for  a  cavity  with  the  parameters:    RQ  =  0.5  ft,  Poo  =  2200  lb/ft    ,  and 

3 
Pr~    1.94  slug/ft    .     Figure  7  is  the  solution  to  the  problem  with  the 

2 
following  initial  conditions:      RQ  =  1.0  ft,  Poo  =2500  lb/ft     ,  and 

fe=  1.60  slug/ft3  . 


B.         ISENTROPIC  CAVITY  OSCILLATIONS 

The  second  problem,  which  assumed  the  cavity  was  filled  with  gas 
and  followed  a  series  of  isentropic  compressions  and  expansions,  was 
approached  in  a   manner   similar  to  the   first  problem.    Starting   with 
R  (  0)   =  R0  and   P  (R,  0)   =   P0  ,    it  was  known  that  for  isentropic 
compression 


P.     ,  ..if 

(5) 


P-  /         \  A 

which   when  applied  to  a  homogeneous  gas-filled  sphere  yields 

This  relationship  was  derived  by  assuming  that  R  <<  a,    where   a   was   the 
speed  of  sound  in  the  gas.    Combining  equations  (3)  and  (6)  it  was  found 
that 


a 


f 


R2JR 


(7) 
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A  dimcnsionless  time  was  defined  as 


{    -        lJ  3    T 

rT~  (8) 


and  combined  with  equation  (7).    The  final  expression  for  cavity  motion 

r 


/ 


JX 


r 


J 


x*[>C(<+i)-<xa'ir    (9) 


where    o<  =     PG/(  V    -1)  poo   and  X  =   R0/R.      Equation  (9)  was  solved 
through  the  use  of  a  numerical  integration  computer  program.    The  cavity 
radius,  represented  by  X  =  R0/R,  was  found  to  oscillate  about  the  equi- 
librium position  X  =  1.      This  response  was  expected,  since  the  cavity 
started  at  an  initial  radius  RQ  ,  and  there  was  no  damping    term  included 
in  the  analysis.    This  departure  from  the  true  cavity  response  was  brought 
about  by  the  assumption  of  isentropic  compression  and  expansion  of  the 
air  inside  the  cavity  and  the  neglect  of  compressibility  effects  in  the 
fluid  outside  the  cavity.     However,  valuable  information  can  be  gained 
concerning  the  frequency  and  peak  pressures  associated  with  the  cavity's 
repeated  expansion  and  collapse.    To  aid  in  the  study  of  these  parameters, 
the  computer  output  was  designed  to  reflect  the  cavity  pressure  fluctua- 
tions with  time.    This  was  accomplished  by  plotting  the  output  in  the 
following  form:    cavity  pressure/ambient  fluid  pressure  versus  time.    This 
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pressure  is  directly  related  to  the  cavity  radius  as  shown  by  equation 
(6).    An  example  of  the  output,  which  was  computed  for  an  initial  pres- 
sure ratio  of  6.67,  is  shown  in  Fig.   8.    This  plot  illustrates  the  cavity 
pressure  oscillation  with  time,  and  further  it  predicts  that  the  peak 
pressure  occurs  when  the  cavity  is  at  its  minimum  radius.     This  fact 
has  been  shown  experimentally  by  Cole    (Ref.    1). 
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Fig.   8       Hydraulic  Ram  Cavity  Pressure  Oscillation 
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IV.     APPARATUS 

The  test  apparatus  consisted  of  projectile  accelerators,  test  speci- 
mens,    a  test  tank,  and  associated  instrumentation  for  evaluating  the 
effect  of  high-velocity  impact  into  liquid  filled  containers. 

A.  PROJECTILE  ACCELERATORS  AND  PROJECTILES 

A  0.22  caliber  rifle  was  used  to  accelerate  the  projectiles.    The 

rifle  was  mounted  on  a  stand,  as  shown  in  Fig.  9,  and  located  about  6.9 

meters  from  the  target  tank.    Impact  tests  were  made  with  projectile 

i 
velocities  in  the  order  of  0.381  km/sec  (1250  ft/sec). 

The  cartridges  were  factory  loaded  0.2  2  caliber  long  rifle  ammuni- 
tion with  a  solid  lead  projectile  weighing  1.95  grams. 

B.  TEST  TANK 

The  test  tank  used  for  the  impact  tests  was  fabricated  from  trans- 
parent plexiglas,  Type  G,   25.4  millimeters  (1  inch)  thick.    The  tank  was 
square  in  shape,  with  the  plexiglas  making  up  the  top,  bottom,  and  two 
sides.    The  front  and  rear   of  the  tank  had  removable  and  replaceable 
walls  into  which  impact  was  made.    A  photograph  of  the  tank  with  its 
support  table  is  shown  in  Fig.    10.     In  this  investigation,  the  removable 
walls  were  made  from  7075-T6  aluminum  sheet  1.27  millimeter  (0.050) 
thick,   secured  to  the  transparent  walls  by  a  clamping  frame. 
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Fig.   10      Test  Tank  Installation 
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C.  RANGE  LAYOUT 

The  basic  elements  of  the  ballistic  range  are  shown  in  schematic 
form  in  Fig.    11.    The  down-range  photograph  (rig.    12)  shows  the  velocity 
sensors  and  the  three  shadowgraph  stations  located  along  the  flight  path. 
These  shadowgraph  stations  were  used  in  conjunction  with  the  velocity 
sensors  to  accurately  determine  the  velocity  and  attitude  of  the  projectile 
passing  down  the  range. 

D.  INSTRUMENTATION 

The  projectile  was  a  chronograph  screen  with  a  5- volt  d.c.  potential 
across  it.    When  the  screen  was  broken  by  a  projectile,  a  pulse  was  sent 
through  a  signal  conditioner,  then  to  the  delay  unit,  and  finally  to  a 
counter.    Two  Monsanto  10 IB  counters  were  used  for  velocity  measure- 
ment.   The  delay  units  were  used  to  trigger  the  shadowgraph  spark 
sources  at  the  proper  time.     Figure  13  shows  the  circuitry  of  the  two 
delay  units  used  in  timing  the  dual  shadowgraph  system  of  the  test  tank. 
The  major  component  of  each  delay  unit  was  a  Signetics  N74123  retrigger- 
able  monostable  multivibrator.     These  units  generate  variable  time  delays 
from  200^cs  to  1700^«s  and  from  Zms  to  18^<*s  ,  respectively.    A  spark 
source  triggering  capacitor  was  built  into  each  delay  unit.    The  signal 
conditioner,   shown  in  Fig.   13,  was  needed  to  compensate  for  the  wide 
range  of  resistances  found  in  the  chronograph  screens.    The   signal 
conditioner  was  basically  a  transistor  switch  using  a  2N3404  transistor. 
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Fig.    13      Delay  Unit  Circuitry, 
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E.         SHADOWGRAPH  SYSTEM 

The  shadowgraph  system  associated  with  the  test  tank  consisted 
of  two  similar  stations.     Figure  14  is  a  schematic  representation  of  the 
apparatus  showing  a  down-range  view.     Each  shadowgraph  station  was 
composed  of  a  spark  source,  a  16-inch  collimating  mirror,  and  a  4-x  5- 
inch  photographic  plate.    The  spark  source  was  located  at  the  focal  point 
of  the  mirror  to  provide  the  necessary  collimated  light.    Polaroid  Type  57 
high-speed  film  was  used  to  record  the  shadowgraphs  through  the  water 
filled  test  tank.    The  complete  shadowgraph  system  is  shown  in  Figs.  15 
and  16. 
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Station  1 


Fig.   14     Schematic  of  Shadowgraph  Apparatus 
(Down-  range  View) . 
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Fig.    16     Dual  Shadowgraph  System 

(Plane  Mirror  and  Film  Holders) 
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V.      RESULTS  AND  DISCUSSION 

This  research  was  conducted  as  a  preliminary  investigation  of  the 
mechanisms  involved  in  hydraulic  ram.     Primarily,  the  phenomena  studied 
were  the  characteristics  of  the  shock  wave  pressure  pulse  generated  in 
water  from  impact  by  a  projectile  and  the  rate  of  energy  loss  of  the  pro- 
jectile in  the  water.     By  the  characteristics  of  the  shock  wave  pressure 
pulse    one    means    the    shape   and   progress    of  the   shock  wave  with  time 
and  the  magnitude  of  the  pressures  at  the  wave  front. 

A  .  SHAPE  OF  THE  SHOCK  WAVE  FRONT 

The  investigation  of  Ref .   3  indicated  that  the  shape  of  the  pressure 
wave  front  generated  in  the  water  by  a  projectile  impact  was  hemispheri- 
cal.    However,  the  data  of  this  reference  were  primarily  limited  to  impacts 
with  spherical  projectiles.     Since  the  typical  threat  to  an  aircraft  fuel 
tank  has  the  standard  bullet  shape,  the  projectiles  used  in  the  investiga- 
tion were  solid  lead  bullets.     The  influence  of  lead  bullets,  which  deform 
as  they  enter  the  tank  and  traverse  the  fluid,  was  explored  and  will  be 
discussed  in  a  later  section  of  this  report. 

The  shape  and  propagation  of  the  pressure  wave  front  resulting  from 
the  projectile  impact  were  obtained   from  the  dual  high-speed  shadowgraphs. 
A  typical  dual  shadowgraph  sequence  of  the  pressure  wave  front  progress 
is  shown  in  Tig.    17.    The  figure  illustrates  the  results  of  the  impact  by  a 
1.95  gram  lead  bullet  at  a  velocity  of  0.374  kilometers  per  second 
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(1226  ft/sec)  into  water  in  the  transparent  plexiglas  test  tank.     For  this 
test  the  projectile  passed  through  a  0.5-inch  diameter,  "pre-punched" 
hole  in  the  entry  wall.    A  strip  of  masking  tape  was  used  to  cover  the 
hole  and  contain  the  water.    A  hemispherical  wave  front  was  produced, 
and  the  center  of  the  hemisphere  remained  at  the  point  of  impact  as  the 
wave  front  expanded.     The  photographs  of  Fig.    17  show  that  there  was 
no  significant  effect  of  the  forward  momentum  of  the  projectile  to  move 
the  wave  front  in  the  direction  of  flight. 

A  commonly  used  model  for  analytical  prediction  of  the  hemi- 
spherical wave  strength  and  motion  is  the  blast  wave  theory.     Blast 
wave  theory  assumes  sudden  release  of  energy  centered  at  the  projectile 
penetration  point.    The  behavior  of  the  fluid  properties  as  a  function  of 
radius  will  not  be  consistent  with  that  model.     See,  for  example, 
Chapter  4  of  Sedov  (Ref.  11). 

An  alternate  and  more  accurate  model  would  be  to  assume  an 
expanding  hemispherical  piston  due  to  projectile  penetration.    The  rate 
of  change  of  piston  radius       \~p      would  be 

Alt    [3     jyi^3 

n   ~Lzir   uj  (10) 

The  quantity    tfv/^i      equals   the  rate  at  which  the   projectile   dis- 
places  fluid    volume  during   penetration. 
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Fig.   17      ShockWaves  Resulting  from  a 
Pre-punched  Wall  Impact. 
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Tests  were  also  performed  with  a  solid  front  wall  on  the  test  tank. 
The  entry  walls  used  were  7075 -T6  aluminum  sheet  1.27  millimeter 
(0.050  in.)  thick.    An  example  of  the  resulting  wave  front  due  to  solid 
wall  penetration  is  shown  in  Tig.    18.    This  set  of  photographs  illustrates 
the  obscuration  of  the  interior  of  the  shock  wave  due  to  the  projectile 
penetration  of  the  entry  wall.     The  shadowgraphs  show  that  the  hemispher- 
ical wave  front  was  produced  as  in  the  "pre-punched  wall  case."      The 
other  wave  front  in  the  water,  which  appears  as  a  straight  line  oblique 
to  the  wall,  was  caused  by  the  propagation  of  shear  and  dilatation  waves 
(Ref .    12)  moving  through  the  metal  at  a  faster  rate  than  the  waves  through 
the  water.    The  rate  of  propagation  of  this  oblique  wave  was  found  to  be 
proportional  to  the  bulk  sonic  velocity  (dilatation  wave)  for  aluminum 
(6  .  2  km/sec) . 

B.  PROGRESS  OF  THE  SHOCK  WAVE  WITH  TIME 

Impacts  with  projectiles  of  various  kinetic  energies  were  made  to 
determine  the  influence  of  kinetic  energy  on  shock  wave  velocity.    The 
average  velocities  of  the  shock  waves  were  determined  from  shadowgraphs 
similar  to  those  in  Figs.    17  and  18.    The  shadowgraphs  were  used  to  ob- 
tain time  and  position  information  for  the  shock  wave  as  it  expanded  and 
propagated  through  the  fluid.     The  slopes  of  the  resulting  time  -  distance 
plots  supplied  the  average  propagation  velocities  of  the  shock  waves. 

Figure  19  is  a  plot  of  the  resulting  shock  front  velocities  for  an 
impact  kinetic  energy  of  135.2  joules.     This    plot   shows   that  the  shock 
front   velocity  decays   to  the   acoustic   speed   in   water  by  a   distance 
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Fig.   18      Shock  Waves  Resulting  from  a 
Solid  Wall  Impact. 
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Fig.    19     Shock  Front  Velocity  Generated  in  Water 
by  Projectile   Impact. 
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of  about  11  centimeters  (4.3  in.)  for  this  impact  energy.    The  acoustic 
speed  for  water  used  in  this  investigation  was  1.5  kilometers  per  second 
(4920  ft/sec). 

C.         PRESSURES  GENERATED  IN  THE  FLUID 

The  pressure  generated  in  the  fluid  by  the  projectile  was  divided 
into  two  regions.    The  first  region  consisted  of  the  pressure  associated, 
with  the  shock  front  propagation  and  was  considered  as  the  shock  phase 
of  the  hydraulic  ram  phenomenon.    The  second  pressure  region  was  con- 
centrated in  the  immediate  vicinity  of  the  projectile  itself. 

1 .         Shock  Front  Pressure 

Pressures  generated  in  the  water  at  the  shock  front  correspond 
to  the  measured  shock  front  velocities  and  were  determined  from  data  given 
in  Table  2.2  in  Ref.    1.     Figure  20  shows  the  pressures  generated  at  the 
shock  front  corresponding  to  the  data  given  in  Fig.    19  of  this  report.    The 
result   shows  that  extremely  high  pressures  are  generated  in  the  water 
near  the  impact  point.     For  example,  impact  with  a  1.9  5  gram  bullet  at  an 
impact  velocity  of  0.372  kilometers  per  second  (1222  ft/sec)  resulted  in 
a  pressure  of         3750         kilograms  per  square  centimeter  (   53,335         lb/sq 
in.)  at  a  distance  of  2.24  centimeters  (0.88  in.)  from  the  impact  point. 
These  pressures  decay  rapidly  due  to  geometric  expansion  of  the  shock 
front  and,  as  in  the  case  of  the  shock  velocity,  reach  ambient  values  in 
a  short  distance . 
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Fig.   20     Shock  Front  Pressure  Generated  in 
Water  by  Projectile  Impact. 
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2 .         Pressure  Field  Around  the  Projectile 

As  with  any  object  moving  through  a  fluid,  the  projectiles 
created  a  dynamic  pressure  in  the  water  due  to  the  stagnation  conditions 
on  their  front  surface.     Dynamic  pressure  is  related  to  the  density  of  the 
fluid  and  the  square  of  projectile  velocity.     Tor  the  velocities  of  0.372 

kilometers  per  second  used  in  this  research,  pressures  of  710  kilograms 

4  2 

per  square  centimeter  (1.01  x  10     lb/in    )  were  found  to  exist  in  the  vicin- 
ity   of  the  projectile . 

Further  insight  into  the  pressure  field  around  the  projectile  was 
gained  by  examining  the  light  diffraction  patterns  illustrated  in  Fig.  21. 
The  concentration  of  light  at  the  front  surface  of  the  projectile  was  caused 
by  a  change  in  the  refractive  index  of  the  fluid  due  to  the  large  increase 
of  pressure  and  density  in  this  region.     Light  is  refracted  toward  the  larger 
refractive  index  which  accompanies  higher  density  and  pressure  in  the  fluid 
A  qualitative  analysis  of  the  entire  pressure  field  around  the  projectile  can 
be  made  by  studying  the  amount  and  direction  of  deflection  of  the  grid  lines 
in  this  area . 

Figure  21  also  shows  the  air  filled  conical  cavity  behind  the 
projectile.    This  cavity  is  a  result  of  the  dynamic  pressure  produced  in  the 
fluid  by  the  projectile. 

D.         PROJECTILE  ENERGY  DECAY  AFTER  IMPACT 

Measurements  were  made  of  the  progress  of  the  projectile  through  the 
water  after  impact.     From  these  measurements  the  projectile  velocity  and 
energy  loss  were  determined. 
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One  set  of  the  shadowgraphs  from  which  the  projectile  velocities 
were  calculated  is  shown  in  Fig.   22.      These  experimentally  derived  veloci- 
ties   and      energies  are  plotted  as  a  fraction  of  their  impact  magnitude 
versus  distance  into  the  tank,  as  shown  in  Tigs.   23  and  24. 

A  detailed  analysis  of  the  velocity  decay  or  energy  loss  of  a  projec- 
tile traversing  a  fluid  would  be  extremely  difficult  because  of  the  many 
variable  parameters.    Assumptions  must  be  made  to  account  for  the  changes 
in  the  projectile's  shape,  drag  coefficient,  and  material  properties.    Because 
of  these  difficulties,  a  simplified  analysis  of  the  projectile  velocity  and 
energy  loss  was  made  and  compared  with  the  experimental  data. 

The  following  drag  equation  was  used: 

Rearranging  terms ,  knowing    V  =   dx/dt,    gives 

-  iV  =    f'UAJx  (12) 

V  z  m 

Assuming  a  constant  deformed  projectile  area  of  0.5  square  centimeters, 
constant  drag  coefficient  of  1.0,  and  constant  projectile  mass,  integration 

yields:  _  ffCpAx 

y  zm 

"       ^  (13) 


V 


Squaring  both  sides  of  this  relation  gives  the  ratio  of  the  projectile  energy 
after  impact  to  the  projectile  impact  energy.    These  calculated  ratios  for 
the  projectile  velocity  and  energy  are  compared  with  the  experimental  data 
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Fig.  22       Projectile  Progress  into  the  Tank. 
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Fig.   23       Projectile  Velocity  Decay  Ratio  as  a 
Function  of  Distance  into  Tank. 
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Fig.   24      Experimental  and  Calculated  Ratios 
of  Projectile  Energy. 
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in  Figs.  23  and  24.  This  comparison  indicates  that  in  spite  of  the  simpli-. 
fying  assumptions,  the  agreement  between  the  calculated  and  experimental 
values  of  the  velocity  and  energy  after  impact  was  fairly  good. 

In  the  experimental  tests,  the  actual  projectiles  underwent  a  con- 
tinual deformation  during  deceleration  in  the   tank.    An  increase  in  frontal 
area  accompanied  the  deformation,  resulting  in  an  average  final  projectile 
area  of  0.5  square  centimeters.    This  was  the  area  used  in  the  calculated 
approximation  of  the  velocity  and  energy  decay  and  constitutes  the  major 
cause  for  error. 

Energy  decay  of  the  projectile  can  be  directly  related  to  the  amount 
of  energy  being  transferred  to  the  fluid  by  the  projectile.     More  importantly, 
the  slope  of  the  projectile  energy  decay  curve  gives  the  rate  of  energy 
transfer  to  the  fluid,  assuming  no  projectile  deformation.     This  energy 
transfer  rate  has  been  established  as  a  major  factor  in  determining  the 
severity  of  hydraulic  ram  (Ref.  7).     The  average  energy  transfer  rate  for  the 
projectiles  tested  was  20.6  joules  per  centimeter  of  fluid  traversed. 

E.  PROJECTILE  DEFORMATION  AFTER  IMPACT  IN  FLUID 

The  projectiles  were  recovered  after  each  experimental  test  and  were 
found  to  display  a  unique  pattern  of  deformation.    The  deformation  phenome- 
non   was  very  consistent  for  the  projectiles  impacting  the  fluid  through  the 
pre-punched  plates.    Three  examples  of  projectiles  after  this  type  of  impact 
are  shown  in  Fig.   25.    The  extent  of  deformation  was  found  to  vary  directly 
with  the  impact  velocity;  however,  all  of  the  projectiles  exhibited   an  iden- 
tical pattern  of  having  one  flat  side  and  one  undeformed  side.    This  pattern 
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Fig.  25      Projectiles  After  Impact  Tnrough 
Pre-punched  Entry  Wall 
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of  deformation  suggests  a  positive  dynamic  stability  for  the  projectiles 
as  they  traverse  the  fluid  after  impact.     Once  initially  flattened,    the  pro- 
jectile maintains  a  fixed  attitude  due  to  its  stability,  and  all  further 
deformation  takes  place  on  the  same  flattened  side.    As  shown  in  Fig.  26A, 
the  dynamic  stability  of  the  projectile  is  maintained  through  an  equilibrium 
condition  between  the  pressure  force,  located  at  the  center  of  pressure, 
and  the  inertia  force,  acting  through  the  center  of  gravity  of  the  projectile. 
Figure  26B  demonstrates  the  effect  of  displacing  the  projectile  from  the  equilib- 
rium    position.    The  center  of  pressure  moves  upward  on  the  slightly  curved 
front  surface  of  the  projectile  with  the  stagnation  point;  this  creates   a 
restoring  moment  which  returns  the  projectile  to  the  stable  position. 

The  projectiles  which  penetrated  a  solid  entry  wall  before  traversing 
the  fluid  did  not  demonstrate  the  previously  discussed  deformation  phenome- 
non.      Two  typical  projectiles  of  this  type  are  shown  in  Fig.    2  7  .    The 
small  aluminum  cap  located  above  each  projectile  in  the  photograph  is  the 
piece  of  the  entry  wall  punched  out  by  the  projectile  as  it  passed  through 
the  aluminum  plate.    A  photograph  of  an  entry  wall  after  a  test  shot  is 
shown  as  Fig.   28 . 

F.         SHOCK  AND  CAVITY  PFIASE  SEPARATION  POINT 

When  the  projectile  first  penetrates  the  wall,  a  strong  shock  wave 
pressure  pulse  is  created.    As  shown  in  Fig.   29,  this  pressure  decays  in 
strength  as'  it  progresses  outward.    At  time   t^   in  Fig.  29  the  pressure 
behind  the  shock  greatly  exceeds  the  stagnation  pressure  on  the  projectile 
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Fig.   26      Dynamic  Stability  of  a  Deformed 
Projectile  in  Water  After  Impact. 
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Fig.  27      Projectiles  After  Impact  Through 
Solid  Entry  Wall. 
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Fig.  28      Solid  Entry  Wall  After  Impact, 
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Fig.   29    Evolution  of  Pressure  in  the  Transition  from  Shock  to  Cavity  Phase, 
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nose.    Comparing  pressures  at    tA  and   t„  ,    it  is  apparent  that  both  the 
pressure  behind  the  shock  and  at  the  projectile  forward  stagnation  point 
are  decaying;  however,  the  decay  of  the  pressure  behind  the  shock  is  much 
more  rapid.    At  time  tc   the  stagnation  pressure  exceeds  the  pressure  behind 
the  shock  wave.    Therefore,  at  some  point  between  t  B  and  tc    the  two 
pressures,  shock  and  stagnation,  are  equal.     One  could  define  this   as 
a  characteristic  time  separating  the  shock  and  cavity  phases  of  hydraulic 
ram.    The  characteristic  time  for  the  impact  energy  of  135.2   joules  shown 
in  Fig.   29  is  32  microseconds. 

A  comparison  can  be  made  to  further  understand  the  differences  between 
the  shock  wave  and  projectile  pressures.    The  stagnation  pressure  created 
on  the  nose  of  the  projectile  as  it  passes  through  the  fluid  can  be  modeled 
by  steady  flow  and  the  equation   (14)  from  Ref.    13, 

Po    .  " 


P 


/  +  ^L(V- 


y-i 

(14) 


while  the  discontinuity  of  the  shock  wave  requires  the  use  of  the  following 
unsteady  flow  equation  (Ref.    13). 

2X         \l/l 


\4.=       (P*;-f)  V+l 


Y        /   p  +Xz±  \  (is) 

For  the  same  fluid  velocity   (V/  c0  -  0.5)    the  pressure  increase  is  19  per 
cent  for  steady  flow  and  84  per  cent  for  unsteady  flow.      This  comparison 


65 


explains   the   large   initial  pressures   of  the   shock   wave   as   well  as   their 
rapid   decay  since  unsteady  flow  is   transitory. 
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VI.     ALTERNATE  APPROACHES  TO  THE  PROBLEM 

As  a  result  of  the  work  reported  here,  the  following  apparatus  and 
test  procedure  are  recommended  for  use  in  the  future  study  of  the  hydraulic 
ram  phenomenon. 

A.  HIGH  SPEED  FRAMING  CAMERA 

The  use  of  a  high  speed  framing  camera  would  be  the  optimum  way  to 
visually  evaluate  the  hydraulic  ram  phenomenon.    A  camera  capable   of 
recording  200,000  frames  per  second  would  be  ideal  for  further  study  of 
the  short-lived  events  which  occur.    This  rate,  a  frame  every  5  microseconds, 
would  supply  a  complete  time-history  for  the  motion  of  the  shock  front  and 
projectile  through  the  fluid. 

B.  PRESSURE  TRANSDUCER 

Pressure  transducers  are  essential  for  the  complete  analysis  of  the 
complex  hydraulic  ram  pressure  field  and  are  recommended  for  use  in  future 
work.    The  pressure  distribution  is  extremely  complex  due  to  tank  wall 
reflections  and  the  continued  expansions  and  collapses  of  the  cavity. 
Recording  apparatus  for  these  pressure  phenomena  require  a  high  frequency 
response  (about  20  kHz  for  cavity  phase  and  1  MHz  for  shock  phase)  and 
are  very  sensitive  to  location  (Ref.  6).     However,  it  is  important  that  this 
type  of  direct  pressure  readout  be  compared  with  previously  obtained 
information. 
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C.  STRAIN  GAUGES 

After  projectile  penetration,  a  bending  moment  is  produced  in  the 
entry  wall  and  propagates  with  the  expanding  shock  wave  pressure  pulse. 
Strain  gauges  would  be  valuable  in  determining  both  the  magnitude  and 
propagation  speed  of  the  bending  moment  and  its  associated  pressure  wave 
front . 

D.  HYDROGEN  BUBBLE  DISPLACEMENT 

In  the  cavity  phase  of  hydraulic  ram,  knowledge  of  fluid  motion 
within  the  tank  is  important.     Hydrogen  bubble  flow  visualization  techniques 
are  well  suited  for  investigations  requiring  determination  of  the  velocity 
field  over  an  extended  region  of  space  in  a  time-dependent  flow.    This 
technique  could  be  used  with  the  transparent  test  tank  and  high  speed 
framing  camera  to  determine  the  pressure  distribution  surrounding  the  cavity 
as  it  expands  and  collapses. 

E.  HOLOGRAPHIC  INTERFEROMETRY 

In  order  to  study  a  material  under  shock  loading,  it  is  necessary  to 
measure  both  the  shock  transit  time  and  the  profile  of  the  transmitted  stress 
wave.    Wall  displacement  information  of  this  type  can  be  obtained  most 
efficiently  through  the  use  of  holographic  interferometry .    The  shocked 
specimen's  free  surface  motion  measured  by  this  technique  could  be  studied 
and  related  to  overall  fuel  tank  response  to  hydraulic  ram. 
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VII.     CONCLUSION 

The  following  conclusions  were  obtained  from  the  investigation  of 
the  shock  wave  and  cavitation  effects  of  the  hydraulic  ram  phenomenon: 

A  ballistic  range  and  a  simulated  fuel  tank  test  apparatus  have  been 
designed,  built,  and  made  successfully  operational. 

The  nonlinear  equations  describing  cavity  oscillation  and  collapse 
have  been  programmed  and  numerical  results  obtained. 

The  rapid  energy  transfer  of  the  shock  phase  produces  piston-like 
energy  release  resulting  in  a  hemispherically  shaped  pressure  wave  propa- 
gating from  the  point  of  projectile  impact.    The  hemispherical  shape  of  the 
shock  wave  was  found  to  exist  for  impacts  through  thin  solid  metal  and 
pre-punched  tank  entry  walls . 

Very  high  pressures  were  produced  in  the  fluid  as  a  result  of  the  pro- 
jectile impact.    A  pressure  of        3750         kilograms  per  square  centimeter 

2 
(     53,335        lb/in    )  at  a  distance  of  2  . 2  centimeters  from  the  impact  point 

was  produced  by  a  1.95  gram  projectile  with  an  impact  velocity  of  0.375 

kilometers  per  second. 

Good  agreement  was  found  between  simplified  analytical  predictions 
and  experimental  data  for  energy  lost  by  the  projectile  as  it  decelerated  in 
the  fluid  . 

Large  dynamic  pressures  were  produced  on  the  projectile  as  it  passed 
through  the  fluid.    These  pressures,  in  the  order  of  700  kilograms  per  square 
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4  2 

centimeter    (10     lb/in    ),    caused    substantial  deformation  of  the    solid 

lead  projectiles  . 

A  characteristic  time  was  defined  which  establishes  the  point  of 

separation  between  the  shock  and  cavity  phases  of  hydraulic  ram.     For  an 

impact  energy  of  135.2  joules  the  characteristic  time  was  found  to  be  32 

microseconds  . 
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APPENDIX  A 

DERIVATION  OF  THE  EQUATION  OF  MOTION 
FOR  THE  HYDRAULIC  RAM  CAVITY 

From  continuity  one  obtains 


which  integrates  to  the  form 


(rzu)  =  0 


a. 


rx    =    F(t) 


Since  there  is  potential  flow 

u   -   ~ 
c)  r 

Integration  yields 

<p(rj)  =  - 

"r«"    .   6  (t) 

Differentiating  with  respect  to  time 

gives 

^4>          ftfi' 

For  an  incompressible  fluid  the  momentum  equation  in  terms  of  potential 
functions  is 


\±tL  ,  L(^\l=cit 


it    r      2  ijf 
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Considering  boundary  conditions  at  infinity    cJ    Q?      /J  t    =       ( 

p 
oo 
and     C  (t)   =    ■■■■         .  Combining  equations  leads  to 

-£-  -  Ik.  +  ftR'^Rft'  _    R* R* 


For  conditions  inside  the  cavity,    r  =  R  ,    and   the   preceding   equation 
becomes 

3-7  P(P.      i)     -      Poo 

AR    *  x  A x   =    ' 


which   is   the   equation   of  motion   for  a    cavity. 
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COMPUTER  PROGRAM 

//FCLM  JCB  (C738,1028,AR22)  ,  »HCLf" 

//  EXEC  FCRTCLGP, REGION. GC=100K 

//FCRT.SYSIN  CD  * 

C 

C 

C      CAVITY  PROBLEM  NUFEER  1 

C 

c 

Q  a*  aj  #  A  5}  3^c  s.'c  Jjt  A  ;£  *  *  A  :£  *  A  $  3fc  5J!  A  *  ^  #  Sj<  *  -  V  #  *  A  -"*  £ $  #  A  ^  *  afe  "V  J)c  V;  &  A  >'<  :}  >!<  *  A  *  A  it  ^  A  A  A  A  A 

c  *  '* 

C  *  * 

C  *  TFIS  PRCBLEM  USES  ThE  NUMERICAL  IMEGRAT1CN  SLB-      * 

C  *  RCUTINE  QATR  TO  SOLVE  THE  NON-LINEAR  CAVITY  ECUATICN  * 

C  *  OF  KjTICN.  PRC6LEP  NUMBER  1  IS  CESIGNED  TC  CONFUTE    * 

C  *  TFE  COLLAPSE  TIME  CF  A  SPHERICAL  VACUUM  FILLEC  HYC-   * 

C  *  RALLIC  RAM  CAVITY. 

C  *  * 

C  *  * 

Q  3)  afr  A  $  ^  $  £  *  *  ^c  *  *  **  *  s}t  *s{«  *  :je  A  ^  »  *  *  *  *  *  *  si*  #  *  *  *  *  *  ^  *  *  *  £*  *  *  =4  *  *  *  *  *  *  *  *  *  *  ^  A 

C 

EXTERNAL  FCT 

DINE  LSI CN  T(20G) , XX 12 CO) ,  AUX(20C) 

WRITE  (6,13) 
C 

C      TFE  INITIAL  CONDITIONS  FOR  PRESSURE,  DENSITY,  AND 
C      CAVITY  RADIUS  ARE  DEFINED. 
C 

PINF=25C0. 

RFC  =  1  .6C 

R=1.0 

A=-l . 0/ ( ( 2 .0/3 . 0 ) * ( P I NF/RHO ) ) **C  . 5 
C 

C       INTEGRATION  PARAMETERS  FOR  SLERCUTINE  QATR  ARE 
C      ASSIGNED. 
C 

XL=C.99CC0CC0 

XU=C.99CCC000 

EFS=0.0C1 

N=LCO 

NCIN=40 

CC  1  K=1,N 


CALL  CATR(XL,XU, EPS , NC I N , FCT  ,  Y ,  IEP,AUX) 


E  =  R*A 

T(K  )=-B*Y 

XX(K)=XL 

XL  =  XL-C  .01CCG000 
1.  CCNTINUE 

XX(  l)  =  l  .OCCO' 

XX(  100) =C.COOO 

CC  2  N=l,100 

kFITE  (6,14)  M,XX(F) ,T(M) 
2   CCNTINUE 
C 

C      SCLUTICN  ARRAY  IS  PREPARED  FOR  PLOTTING  USING  SUE- 
C      RCUTINE  CRAW. 


REAL*8  TITLEU2) 
READ  (5,16)  TITLE 
REAL  LABLE/'     '/ 

CALL  DRAW  (  100  ,  T  ,  XX  ,0  ,  0  , LABLE  ,  T  I  TL  E  ,  0  .0  ,  C  .  C  ,  C  ,  1  ,0  ,  C  ,  cj  , 
6,0, LAST) 

13  FCPMAT  (• L% /////, T23, «R/RZERG«  ,T42, 'TIFE  TC  CAVITY  CC 

LLAFSE1 ) 

14  FCRVAT  (/,Tlb ,I3,T23,F8.6,T48,F8  .6  ) 
16  FCRfAT  (6A6) 

STCF 
ENC 
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FLNCTICN 
FCT=X**1 

RETURN 

Er^C 
//CC  .SYSIN  OD 
HGLNfD.   TIME 
TO  A 


FCT(X) 

5/(  L.0-X**3. C)**0 


TO  CAVITY  CCLLAPSE  FPCM  RZERO 
RADIUS  R  (PLCTTED  AS  R/RZERO) 
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//H 

// 
//F 

C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 

c 
c 

c 


CLM    JCB    (0738, 1028, AR22) , «hOLM»  ,  T I K  E  =  4 
EXEC    FCkTCLGP,REGI3N.GC=10  0K 

CRT  .SVSIN    CD    * 


CAVITY  PROBLEM  NUMBER  2 


3  ^  $  $  *  if  *  *  if  it  if  if  if  if  sjt  ^  5j  ^  #  #  jJi  *  $  *  ^  :*:  £  ifc  V;  k.  *.  9  •$:  :£  >'<  *  $  if  $  if  $  it  *  r$  >>  ^  :£  #  A  £  if  *  3{"fr  #  *  ~ 

*       ThIS    PRCBLEM    USES    THE    NUMERICAL     JMEGRATICN    SUB- 

N    * 
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THIS    PRCBLEM    USES    THE    NUMERICAL     IMEGRATICN    SUB- 
RCUTINE    QATR    TO    SCLVE    THE    NCN-LINEAR    CAVITY    ECUATICI 
CF    M3TICN.     PRCBLEM    NUMBER    2     IS    CESIGNEQ    TC    COMPUTE 
CSCILLATCRY  MOTION  OF  A  GAS  FILLEC  HYDRAULIC  RAM 
CAVITY.     THE    OUTPUT     IS    A    PLOT    OF     ThE    CAVITY     PRESSURE       * 
R  A  T  1 0  V  S  T  I M  E  . 

EXTERNAL    FCT,F 

DIMENSION    T(500)  ,P(50G)  ,AUX(500)  ■ 

COMMON  GAMMA, ALPHA 

WRITE  (6,19) 

THE     INITIAL    CONDITIONS    FOR    PRESSURE    AND    FRCFERTIES 
CF  THE  GAS  ARE  DEFINED. 

GAMMA=4  .C/3.C 
PP£TI0=6.667 

•TIME'     IS    THE    CUMULATIVE    TIME    CCLMER. 

TIME=0.C 
ITMAX=5C 

N     IS    THE    SOLUTION    ARRAY    COUNTER. 

N=l 

TEE    FIRST    ROOT    OF    TFE    CAVITY    EQUATICN    OF     FCTICN     IS 
FCLKD    USING    SUBROUTINE    ZFALSE. 

>L  =  C30 
XR=C.50 

ALFHA=PRATIC/(GAMMA-1  .0) 

CALL  ZFALSE  .  ( F , 0 . CO  1 , 3 , XL , XR , XA P F  ,  ITMAX ,  I  E R  ) 

FF=PRATIC 

PRCCT=XAPP**(3.0*GAMMA) 
PPATIO=PRCGT*PRAT IC 

INTEGRATION  PARAMETERS  FOR  SLEROLTINE  QATF  ^RE 
ASSIGNED. 

XLCW=0.999 

XL  =  C999 

DC    2    1=N,5C0 

IF     (  I. EC. 311)    GO    TC     IC 

CALL  GATP(XLCW,XG,C.0C1,40,FCT,Y,IFR,AUX) 

T(I  )=Y  +  TIME 

P(I)=XLCW**(3.0*GAMMA)*FP 

XLCW=XLCW-C.020 

IF  (XLOk.LT.XAPP)  GC  TC  3 

CCNT  INUE 

TIME=T(  I) 

N=I 

WRITE  (6,20)  I,XAPP,T  (I  ),P(  I)  ,PP 
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XL=2.20 

THE  SECOND  ROOT  QF  THE  CAVITY  EOLATION  CF  NCTICN  IS 
FCLND  USING  SLBRGLTINE  ZFALSE. 

XF=2.40 
ALFhA=PRATIO/(GAMMA-l.G) 

C/>LL    ZFALSE     (  F  ,  0  .  0  0  1  f  2  ,  XL  ,  XR  ,  XAPF  ,  1  TMAX  ,  I  EF<  ) 

PF  =  PRAT  IC 

FRCCT=XAPP**(3.0*GAMMA) 

PPATIO=PRCCT*PRATIC 

INTEGRATION  PARAMETERS  FOR  SUBROUTINE  QATP  /IRE 
ASSIGNED. 

XLCW-1.C01 

XL=1.001 

CC  6  I=N,500 

IF  (1. EC. 311)  GO  TC  10 

CALL  QATRJXLQW, XU,C.0C1,40,FCT, Y,IER,AUX) 

T(  1)=Y+TIME 

P(I  >=XU**(3.0*GAMMA)*PP 

XL=XL+C.02C 

IF  (XU.GT.XAPF)  GC  TO  7 

6  CONTINUE 

7  TINE=T(I) 
N  =  I 

WRITE  (6,20)  I,XAPP,T  (I ) ,P( I  )  ,PP 
GC  TO  1 
IC  WRITE  (6,12) 
CC  11  M= 1,310 
WRITE  (6,14)  f",T(F)  ,P(P) 

11  CONTINUE 

SCLUTIGN  ARRAY  IS  PREPARED  FOR  PLOTTING  USING  SUE- 
RCCTINE  CRAW. 

REAL* 8  TITLE (12) 
RE  AC  (5,15)  TITLE 
REAL  LAeLE/'     ■/ 

C^LL  DRAW  (300, T , P , C , 0 , LABL E , TITLE, 0. 0, 0, C,C  iC,  0,0,6, 8 
,0,LAST  ) 

12  FCPMAT  (M1  ,/////, T25, 'NON-DIMENSICNAL  T  I  N L e  ,  T50 ,  •  C AVI 

TY  PRESSURE/PRESSURE  AT  INFINITY') 

14  FCRNAT  (/,T18,I3,T27,F1C.4,T53, F1C.4) 

15  FORMAT  (6A8) 

19  FCPNAT  («  1*  ,///,T22,  •  I  «  ,T31,  «XAFP  •  ,T43,  M  INE«  ,T55,  «F(I 

)  «  ,T68,  •PPM 
2C  FCPNAT  (/,20X,I3,4  (2X,F10.4)  ) 

STCP 

END 


FLNCTICN  F(X) 

CCPNGN  GAMMA, ALPHA 

F=-(ALPhA*X**( 3.0* GAMMA) )+((ALPPA+1.0)*X**2.C)-l.C 

RETURN 

END 

FLNCTICN  FCT(X) 
CCMMCN  GAMNA,/>LPHA 

PCT=(X**(-2.0) )*(-(ALPHA*X**(3.0*GAMMA)  )+  (  (ALPHA+I  .O* 
X**3.0)-1.0)**(-0.5) 

RETURN 
END 
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S16RCUTINE  ZFALSE  (  F  ,  EPS  ,  NS  IG  ,X  L  ,  X  R  ,XAPF  ,  11  f'A  X  ,  I  BP  ) 


C 

c 

c 


1C 


15 


2C 


IEP  =  0 
IC  =  0 

Ef-SP  =  10.** 
FXL  =  F(XL) 
FXR  =  F(XR) 
IF  (FXL*FXF) 

IER  =  129 
GC  TO  4C 

XAPP  =  XR 

IF  (FXL  .EG. 

GG  TO  4C 

CCNPLTE  APPROBATE  PCCT 
XAPP  =  XL  +  FXL-(XR-XL}/ (FXL-FXR) 
FXAPP  =  F(XAPP) 

IPS)  GO  TO  2C 


(tNSIGJ 
15, L0f5 

0.0)  XAPP  =  XL 


TERMINAL  ERROR 
FXL  CF  FXR  =  G 


IF  (ABS(FXAPP) 
GG  TO  40 


,GT 


IF  (FXAPP^FXL  .GT. 

XP  =  XAPP 

FXR  =  FXAPP 

GC  TO  3C 

XL  =  XAPP 

FXL  =  FXAPP 


DETERMINE  wPETFER  APFPGXI 
LIES  BETWEEN  XAFF  AND  XL 
ANC  XP 
CO)  GO  TO  25 


IF 

GC 


(XR-XL 
TO  4C 


DETERMINE  IF 
.GT.  EPSP*ABS(XR)  )  GO  TC  35 


25  IC  =  IC+1 

IF  (IG.LE.50) 

IEP  =  130 
AC     ITNAX  =  IC 

IF  (IER  .NE. 

GC  TO  9005 
90CC  CCNTINUE 

CALL  UERTST  < IER ,• ZFALSE  ■  ) 
9005  RETURN 

ENC 


(XR-XL  .GT. 

H^X  ITERAT 


GO  TC  15 

0)  GC  TO  9000 


SLEROUTINE  UERTST ( I ER, NAME ) 


IC 


DIMENSION 
INTEGER*2 
INTEGER 
EQUIVALENCE 


CATA 


* 


CATA 

IEP2=IEP 
IF  (IER 2 

IER  1  =  4 

GC  TO  2C 
IF  (IER2 

I  ER  1=3 

GC  TO  20 
IF  (IER2 

IEP1=2 


ITYP 


IEIT 

PRINTR 


ITYP (5,4) ,  IBIT(4) 

NAME  (3) 

WARN,WARF  ,TERN,PRINTR 

(IEIT(l),WARN),(iaiT(2),kAPF),(IEIT 
/•WARN*  ,  '  ING  '  r  '      *  »  ' 
«  UARN'  , '  ING(  '  , •WITH1  ,  ' 


',  ' 


•NGN-' ,  «DEF  I  c  ,  'NED  •  ,  « 
/  22,64,  128, C/ 
/  6/ 


FIX'  ,  '  ) 
'  t  « 


GE.  WARN)  GO  TO  5 


LT.  TERN)  GG  TO  10 


LT.  WARF)  GC  TO  15 


NON-CEFINED 


TERMINAL 


WARNINGCWITP  FIX) 
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15 

2C 


GC  TO  2C 
I  ER  1=  1 
IER2=IER2-IBIT( IER1J 


WARNING 
EXTRACT 


N 


PRINT  ERROR  NESSAGE 
WRITE  (PRINTR,25)  ( I T YP ( I , I ER  1 )  ,  I  =  1 ,  5 )  , NA f E  ,  I  ER2 
25   r-CPVAT(«  ***  I  M  S  L(UERTST)  ***   •  ,  5A4  ,  4X  ,  2A2  ,  4X  ,  I  2  ) 
RETURN 
ENC 
//GC  .SYSIN  DO  * 
hCLNfD.  HYDRAULIC  R  MA    CAVITY 
PRESSURE  OSCILLATION 
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